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Abstract 

The correspondence of the Gauss-Bonnet (GB) and its modification (MGB) models of dark 
energy with the standard and generalized Chaplygin gas-scalar field models (SCG and GCG) have 
been studied in a flat universe. The exact solution of potentials and scalar fields, which describe the 
accelerated expansion of the universe, are reconstructed. According to the evolutionary behavior 
of the GB and MGB models, the same form of dynamics of scalar field and potential for different 
SCG and GCG models are derived. By calculating the squared sound speed of the MGB, GB 
model as well as the SCG, GCG, and investigating the GB-Chaplygin gas from the viewpoint of 
linear perturbation theory, we find that the best results which is consistent with the observation, 
may be appeared by considering the MGB-GCG. Also we find out some bounds for parameters. 
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I. INTRODUCTION 


Astrophysical data which is out coming from distant la supernova [l-3|, Large Scale 
Structure (LSS) 4, 5| and Cosmic Microwave Background (CMB) 0.Q...d.catethato„. 
universe undergoes with an accelerating expansion. This kind of expansion may be arisen 
by a mysterious energy component with negative pressure, so called, dark energy (DE). 

However in the last decades, other models based on modified gravity 
{F{R), F{G), F{R, (j), X), F{T),...) have been proposed that have given another de¬ 
scription of acceleration expansion of the universe. In these models, many authors have 
showed that all models of DE can be resolved by modifying the curvature term R (Ricci 
scalar) of Einstein-Hilbert action with another curvature scalars such as any scalar function 
of R, Gauss-Bonnet term (G), torsion (T), scalar-tensor (X,0) and etc. (details are in 
Ref. js-SlI and references there in). Even, some authors found that, the early inflation, 
the intermediate decelerating expansion and late time acceleration expansion, could be 


described together in one model 


22 |. 


Lately, among many models of DE, dynamical models, which are considering a time 
dependent component of energy density and equation of state, have attracted a great deal 
of attention. Also, among many dynamical models, ones that represented by a power 
series of Hubble parameter and its derivative {i.e. F[, F[F[, F[^,...) have been interested 
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24| . Also authors 


m 


25 


29| have shown that terms of the form H^, HH^ and can 


be important for studying of the early universe. Hence, it would not be some thing strange 
to consider a DE density proportional to the Gauss-Bonnet (GB) term which is invariant in 
4-diniensional. Besides, in geometrical meaning, the GB invariant has a valid dimension of 
energy density 301. Also authors in {m, 321 showed that a unification between early time 
inflation and late time acceleration in a viable cosmology can be described by a coupling 
between GB term and a time varying scalar field 331. 


The other successful model of D 
model (SGG), first proposed by 


34 


is Ghaplygin Gas model. The standard Ghaplygin Gas 
, regards as a perfect fluid which plays a dual role 


in the history of the universe: it behaves as dark matter in the first epoch of evolution of 
the universe and as a dark energy at the late time. Unfortunately this model has some 
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So Generalized 


inconsistency with observational data like SNIa, BAG, GMB 
Ghaplygin Gas (GGG) 40] and Modihed Ghaplygin Gas (MGG) models 41-^| have been 
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introduced in order to establish a viable cosmological model. It would be beneficial to 
study any relationship between SCG model and its modification while DE density behaves 
like GB invariant term as mentioned above. In this paper we would show that it leads 
interesting cosmological implications. 

As we would show in this paper, the EoS parameter of GB DE model on its own does not 
give rise to phantom phase of the universe. Besides in [3^, author shows that presence 


of matter drastically converts Friedmann equation into a nonlinear differential equation 
which alters the behavior of the EoS parameter which can lead to Wo ~ —1.17 and allows 
for quintom behavior. However, in this paper, we incorporate GB dark energy density with 
a SGG component without adding any matter content. In addition, corporation GB or 
MGB with different GG models {i.e. SGG and GGG) would be help full in order to obtain 
exact solution for scalar field and potential and would relieve us in order to determine 
some bounds for free parameters of models. So considering the cosmological solution for 
different compositions of GB and GG models could show the importance of each one. Also, 
we would succeed in the frame work where = SttG = = 1 and in the natural unit 

where {h = c = 1). 

The outline of this paper is as follows: In next section, we introduce the GB dark energy 
and calculate the deceleration and EoS parameters. Then, in subsections 2.1 , 2.2 and 2.3 
we investigate corporation GB with SGG and GGG, in turn and then scalar field and scalar 
potential are obtained by exact solution. In section III, the same procedure has done for 
MGB energy density. In section IV, we would investigate Adiabatic Sound Speed, which 
is one of the critical physical quantity in the theory of linear perturbation. In section V, we 
discuss on behavior of scalar field, scalar potential and deceleration parameter versus x for 
GB and MGB models and we gain some bounds for free parameters of models. Finally, we 
summarize our results in Sec. VI. 

II. GAUSS-BONNET DARK ENERGY IN A FLAT UNIVERSE 

The energy density of GB-DE is given by 


Pd = oiQ 


( 1 ) 
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where a is a positive dimensionless parameter 30|]. Gauss-Bonnet invariant Q is topological 
invariant in four dimensions and may lead to some interesting cosmological effects in higher 
dimensional brane-world (for a review, see 4^). It is dehned as 


Q = R^- 


puri'y 


( 2 ) 


where R, R^j.^ and R^j,vm scalar curvature, Ricci curvature tensor and Riemann curvature 
tensor, respectively. In a spatially flat FRW universe 


d^s = —dt^ + (t) sin^ dd0^] 

the Eq. ([T]) takes the form 


Pd = 24aH^ [H^ + H). 


( 3 ) 


( 4 ) 


By using the energy density Pd, without any matter component, the Friedmann equation in 
flat universe in reduced Planck mass unit {SttG = h = c = 1) is 

= ^pD = 8aH^ [h^ + id) . (5) 

Dehning the e-folding x with dehnition x = Ina = —ln{l + z), where z is the redshift 
parameter and using d/d{x) = -^d/d(f), we get the following differential equation 




2 dx 


8a 


which immediately gives the solution 

H{x) = 


(n-ee-2"). 


8a 


( 6 ) 


( 7 ) 


The parameter ^ is a constant of integration which is obtained by ^ = SaH^ — 1. Also it 
gives a = (1 -|- ^)/{SHq). Using the continuity equation 


Pd + 3H {1 + wd) Pd = 0 

and Eqs. (Ilj),(l5]), the equation of state (EoS) parameter yields 

Pd 


Wd = -1 


= -1 


2 H 


= -1 


1 )- 


( 8 ) 


( 9 ) 


?>HpD ' 3dd2 3'8add2 

It is more preferable to write above equation in term of e-folding, x. Hence, by using Eq. 
o. the EoS parameter can be rewritten as 


Wd = -1 + 77 




-2x 


3 Vl + ee 


-2x 


( 10 ) 
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We see that the constant ^ plays a crucial role in the behavior of the EoS parameter. For 
.^ = 0 (he. SaH^ = 1), the EoS parameter for ACDM model {wa = —1) is retrieved. For 
^ > 0 the expanding universe accelerates in quintessence phase (—1 < wd < —1/3). Using 
Eqs. ([S]) and ([7]), the deceleration parameter is calculated as 


H 1 1 


( 11 ) 


Since a and Hq are positive parameters, so ^ always must be greater than —1. Therefore, 
the deceleration parameter is always negative except for — 1 < ^ < 0. In this way, the 
universe which is characterize by GB dark energy model could not exhibit a transition from 
deceleration to acceleration phase for ^ > 0, against what we expect from observations. 


A. Gauss Bonnet Standard Chaplygin Gas 


The SCG is a perfect fluid with an equation of state as 

PscG = -, (12) 

P 

where p, p and A are pressure, energy density and a positive constant respectively. By 
substituting Eq. flT^ into the continuity equation (0), the energy density immediately solved 


PscG = 'JA + Be 


(13) 


where B is an integration constant 
the energy density and pressure are dehned 


45| . Using the standard scalar held DE model in which 


and equating pscG 


P<f> 

Pcj> 


+ r (<#■) 
- r (^) 


-A 

\/A + ’ 


(14) 

(15) 


and pscG = P<i>^ the scalar potential and kinetic energy term of 


SCG model are given as 


U(0) 

= 


2A + Be-^^ 
2y/A + Be-^^ 

B^-dx 

A + Be-^^' 


Also the EOS parameter becomes 


p A 

WSCG = - = 

P 


(16) 

( 17 ) 
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A + Be-^^ 


( 18 ) 














Equating the energy densities (he., pscG = Pd) and EoS parameters (he., wscg = wd), 
after using the Friedmann equation ([5]), constants A and B immediately given by 


A = 


(8a) 


[(2 + ee-'")'-l], 


B = e 


6x 




A 


and hence the scalar potential and kinetic energy term rewritten as 


^ = h/^. 

^ 2 V a 

By inserting (j)' = (f)/H , where prime means derivative with respect to x - 
ential equation fl22]) gives the normalized scalar held (0 = 1 at present, x 
X as 

^ ^ x/2, ( 1 + 2^6-2" + 2Jie-^-(l + ee-2-) 

(h = 1 -In ' 


1 + 2 ^ + + ^) 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 

In a, the differ- 
: 0) in terms of 

(23) 


It is easy to see that from Eq. d?]), at present, we must have 1 + .^ > 0 and from fl23|l . it 
must be required that ^(1 + ^) > 0. Therefore in this model, we must have ^ > 0. As it 
is shown in Fig. [H the normalized scalar held grows up to a saturated value at late time 
in such a way that this value exceeds for larger values of Also Eq. fl?Ill shows that the 
universe goes to a stable equilibrium at inhnity where E(cx)) = 3Hq/(1 + ^) and from (12^ . 
the scalar held reaches to 0(oo) = 1 + (■\/2/2) In + 2,^ + 2yJ^(l + O) • 


B. Gauss Bonnet Generalized Chaplygin Gas 


The equation of state of generalized Chaplygin gas (GCG) dehned as [46 1 

A 


P = 


P' 


,<5-1 ’ 


(24) 


where A is a constant and 1 < 5 < 2. For 5 = 2, it reaches to SGG model. The energy 
density, similar to previous case, is given by 


PGCG= 


(26) 


and the scalar held model gives energy density and pressure of GGG as 

Pt^lp + V{4,) = {A + Be-^‘A, 


(26) 
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I - ^^O.I - ^=0.2 - ^^0.4 



X 


FIG. 1: behavior of normalized scalar field versus e-folding x for some values of ^ > 0 


-V{4') = -A (A + 


S-1 

s 


(27) 


After forward calculation, three quantities: the scalar potential, kinetic term and EoS pa¬ 
rameter are given by 

V(x) = 


2A + Be 


— 


2(A + 5e-3'5^) « 

'2 

(A + Be-3'5^) « 
p A 

WccG = - = - JTjrBpjE- 

Also same as previous, the constants A and B reconstructed as 


A = 


3 + 


-2x 


(8ay 


[3 (1 + 


B = e 


3Sx 


8a 


(1 + ee-^") 


5-1 


A 


and the potential and dynamics of GB-GCG can be written as 

3 + 2^6-2^ 


V x) = 


8a 


1 

2\ a ' 


(28) 

(29) 

(30) 

(31) 

(32) 

(33) 

(34) 


As it is seen, the potential and dynamics of GB-GGG are not a function of parameter 6 and 
are exactly similar to previous case (see Eqs. fl2l|) and ([22])). Therefore the reconstructed 
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scalar potential and scalar field obtained by previous Eqs. (1^ and ([23]). It is worthwhile to 
mention that both models that we have been studied, encourage with an essential problem. 
Despite of observational predictions, the phase transition between deceleration to acceler¬ 
ation expansion did not happen in GB-DE model. Therefore we will study on the MGB 
model, which may alleviate this problem. 


III. MODIFIED GAUSS BONNET DARK ENERGY 


The energy density MGB has been dehned by 


PD = + \H), 


(36) 


where 7 and A are dimensionless constants 
flat universe gives 


30[ |. The Friedmann equation in dark dominated 


dH^\ 




and the Hubble parameter given by 


H{x) = ,/i(l +pe ’I"), 


(36) 


(37) 


where rj is an integration constant which is obtained by 7 = — 1. The EoS parameter 

becomes 


( 27a: 

r/e ^ 


(38) 


and deceleration parameter is obtained as follows 


27a: 

r}e~~ 


H 7 

Q ~ 1 rT2 ~ ^ T I 1 _22£ 

n A \ I rje \ 


(39) 


For positive values of 7 and A, from Eq. flTT]) . it is easy to see that rj must be always greater 
than —1 and from Eq. fl39|) . a transition from deceleration to acceleration is expected 
provided that 7 > 0. Detailed discussion were transferred to section jVj 


A. Modified Gauss Bonnet And SCG 

Now we want to investigate on the correspondence between MGB-SGG models and recon¬ 
struct the potential and dynamics of scalar held. Same as before, equating energy densities 







[i.e., Eqs. 


and ([35])) and EoS parameters (ne., ffTSD and fl38D ). yield 
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A = 

Y 

B = 


/ 2~fx \ 

\l + r]e > j 


^ 2'y ^ 2-fx 

l + (l-^)r/e A 


3 , 2yx , 

-(1 + rie A ) 

7 


A 


By snbstitnting A and B in Eqs. flTB]) and ffT7|l . we find 


V(x) = - 
7 


O' , 2jx 

1 + (1 - ^)rje A 


2?7 _ 27 a; 


(40) 

(41) 


(42) 

(43) 


which immediately gives the normalized scalar field as 


0 = 1 



2'7a; 

1 + 2Tje 5r + 2 


27CC , 27CC 

rje (1 + ije ^ 


1 + 2 p + 2 -^p(l + ?]) 


(44) 


The behavior of scalar field in this model is the similar to GB-DE model as discnssed in Sec. 


B. Modified Gauss Bonnet And GCG 


As previons, the constants A and B are 


A = (-) 1 + r/e 

7 


27 a; 


<5-1 




B = e 


3Sx 


3 _27£ \ A 

— {l + rje A ) — ^ 

7 


and the potential and dynamics of MGB-GCG are given by 


V{<I>) = - 

7 


1 +( 1 -£)V 


27 a: 


2fj 

T' 


(45) 

(46) 


(47) 

(48) 


which are exactly similar to fl42|) and fl43|) in previons model. Therefore the behavior of 
normalized scalar field and potential are the same as MGB-SGG model. 
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IV. ADIABATIC SOUND SPEED 


Investigation of the sqnared of sonnd speed, would help us to determine the growth 


of perturbation in linear theory 


47l |. The sign of v’l plays a crucial role in determining the 


stability of the background evolution. Positive sign of shows the periodic propagating 
mode for a density perturbation and probably represents an stable universe against pertur¬ 
bations. The negative sign of it shows an exponentially growing/decaying mode in density 
perturbation, and can show sounds of instability for a given model. The squared of sound 
speed is dehned as 

2 _dP _P 

dp p 

In a dark dominated flat universe, it can be written as 

H 


(49) 


n" = -1 - 


HH 


(50) 


and it immediately gives a constant squared of sound speed for GB-DE as = —1/3. 
Therefore it may reveal an instability against the density perturbation in GB-DE model. 
For MGB-DE, Eq. (|5n|l gives 


2 1 , 27 

" =-^+3A' 


(51) 


It shows that can be positive provided that y/A > 3/2. Thus an stable DE dominated 
universe may be achieved in this model. In the next section we would improve this bound 
for 7 /A in a proper way. 


V. DISCUSSION 


We are interesting to focus on MGB-DE model. At hrst, we start with Eq. (13^ and 
plot the deceleration parameter with respect to x in Fig. [2l It shows that the deceleration 
parameter transits from deceleration (g > 0 ) to acceleration (g < 0 ) in some point at the 
past. The parameters 7 and 7 /A play a crucial rule for this point. As 7 or 7 /A adopt 
bigger values, the transition point approaches to present time. By choosing the best values 
oLSoJ, rsj — 0.6) and inflection point as {x ~ —0.5) which has been parameterized recently 
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50|, we obtain some bounds for rj and 7 /A as follow 


0 < 7 < 2.5 


1.5 < Y < 3 
A 


(62) 
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- r|=l .5, y/>.=1 .5 - t|=2, y/X= 1 .8 - r|=0.5, 'ilX=2 

- r|=0.8, y/X=1.8 _ 



FIG. 2: The behavior of deceleration parameter q versus e-folding x for various rj and 7 /A. The 
transition from deceleration to acceleration was happened around x ~ 0.5 


- y/X=2.5 - y/X=2 - y />^= 1-5 



FIG. 3: behavior of V{4>) versus e-folding x for various 7 /A and 7 = 1.5 

Using Eq. fj^ for MGB model, we plot V{(j)) = 7 ^( 0 ) versus x for different values of 7 /A 
and 7 = 1.5 in Fig. [3l This figure shows that as time goes, U((/)) is decreasing to small 
values and the potential will reach to a constant at inhnity. In addition, by increasing the 
ratio of 7 /A, the tracking potential adopts bigger values at future. 

As it is seen, the potential describes a tracker solution. According to the quintessential 
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tracker solution, our universe undergoes a phase from r(; = 0 tota = —1 and the effective 


EoS is Weff = —0.75 The huge advantage of the tracker solution is that it allows 

the quintessence model to be insensitive to initial conditions 52|. So we use this feature 
in order to improve obtained bounds of parameters. In this way, Eqs. fl55]) and flT7|l . for 
matter dominated universe {w = 0), leads to E(0) = 8/(27 — 3A). On the other hand the 
quantity V{(j)) for quintessence barrier {w = —1/3) reach to E(0) = 2/(7 — A), so that the 
value 7 /A = 1 is illegal. It is also consistent with what we got from investigation of the 
deceleration parameter. Finally, the potential might give a tracking solution provided that 
1.5 < 7 /A < 3. 


VI. CONCLUSION 

In this paper, the reconstruction of GB-DE and some variety of Chaplygin gas have been 
studied. We obtained exact solutions for reconstructed scalar held and its potential in each 
models (GB-SCH, GB-MCG, MGB-SCG and MGB-MCG). According to cosmological pre¬ 
dictions and historical evolutions, some models should be rejected (he., models combined 
with GB-DE) and another models which have been combined with MGB-DE can be per¬ 
mitted to express the evolution of the universe. The equation of state and deceleration 
parameters for both GB and MGB models were calculated. In GB-DE model, the deceler¬ 
ation parameter was always negative except for — 1 < ^ < 0. This fact was shown that a 
transition from deceleration to acceleration expansion could not have happened in the past 
that is contrary to the facts of cosmology. Also it was easily shown that the EoS parameter 
in GB-DE model would not ever reach to phantom phase (he. wd < —1)- We showed that 
for = 0 (he. SaH^ = 1), the EoS parameter for AGDM model was retrieved. Investigation 
on the squared of sound speed, revealed an instability of model against density perturbation 
in GB-DE model. 

In MGB-DE model, we found that the transition from deceleration to acceleration is 
permitted just for a limited range of values of rj and 7 /A. Ghoosing the best values for 
deceleration parameter at present and deflection point, according to observations, some 
bounds of 0 < 7 < 2.5 and 1.5 < 7 /A < 3 were obtained. We showed that by redehning 
7 U( 0 ) = D( 0 ), the scalar potential decreased to smaller values and will reach to a saturated 
constant at late time. Our investigation on U(0) for two phases, matter dominate and 
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quintessence, showed that 7 /A, could not take two values 1 and 3/2. 

It will be interesting to hnd the constraints of these models against the data of cosmo¬ 
logical observations and structure formation. We hope to discuss these issues in the future. 
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